Extended Dyer-Roeder Approach Improves the Cosmic Concordance Model 
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A new interpretation of the conventional Dyer-Roeder (DR) approach by aUowing hght received from 
distant sources to travel in regions denser than average is proposed. It is argued that the existence 
of a distribution of small and moderate cosmic voids (or "black regions") implies that its matter 
content was redistributed to the homogeneous and clustered matter components with the former 
becoming denser than the cosmic average in the absence of voids. Phenomenologically, this means 
that the DR smoothness parameter (denoted here by as) can be greater than unity, and, therefore, 
all previous analyses constraining it should be rediscussed with a free upper limit. Accordingly, by 
performing a statistical analysis involving 557 type la supernovae (SNe la) from Union2 compilation 
data in a flat ACDM model we obtain for the extended parameter, as ~ 1.261q54 (Icr). The effects 
of ue are also analyzed for generic ACDM models and flat XCDM cosmologies. For both models, 
we find that a value of a_B greater than unity is able to harmonize SNe la and cosmic microwave 
background (CMB) observations thereby alleviating the well known tension between low and high 
redshift data. Finally, a simple toy model based on the existence of cosmic voids is proposed in 
order to justify why ue can be greater than unity as required by Supernovae data. 

PACS numbers: Dark energy, cosmic distance, supernovae, inhomogeneities 



I. INTRODUCTION 

The accelerating cosmic concordance model (flat 
ACDM) is in agreement with all the existing observa- 
tions both at the background and perturbative levels. 
Hovifever, vifhile more data are being gathered, there is an 
accumulating evidence that a more realistic description 
beyond the "precision era" requires a better comprehen- 
sion of systematic effects in order to have the desirable 
accuracy. 

Local inhomogeneities are not only possible sources 
of different systematics, but may also be signalizing for 
an intrinsic incompleteness of the cosmic description. 
This occurs because the Universe is homogeneous and 
isotropic only on large scales (> 100 Mpc). However, on 
smaller scales, a variety of structures involving galaxies, 
clusters, and superclusters of galaxies are observed. Per- 
meating these structures there are also voids or "black 
regions" (as dubbed long ago by Zel'dovich jlj) where 
galaxies are almost or totally absent as recently suggested 
by the N-body Millenium simulations This means 
that statistically uniform cosmologies are only coarse- 
grained representations of what is actually present in the 
real Universe. As a consequence, the description of light 
propagation by taking into account such richness of struc- 
tures is a challenging task to improve the cosmic concor- 
dance model, but the correct method still remains far 
from a consensus [1, i, S S 0, il • 

Zel'dovich [p, followed by Bertotti Gunn 11 1 and 
Kantowski [l^l were the first to investigate the influence 



of small-scale inhomogeneities in the light propagation 
from distant sources. Later on. Dyer and Roeder (DR) 
assumed explicitly that only a fraction of the average 
matter density must affect the light propagation in the 
intergalactic medium. Phenomenologically, the unknown 
physical conditions along the path, associated with the 
dumpiness effects, were described by the smoothness pa- 
rameter: 
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where ph and pd are the fractions of homogeneous and 
clumped densities, respectively. This parameter quan- 
tifies the fraction of homogeneously distributed matter 
within a given light cone. For a — (empty beam), 
all matter is clumped while for a = 1 the fully homo- 
geneous case is recovered, and for a partial dumpiness 
the smoothness parameter is restricted over the interval 
[0, 1]. The reader should keep in mind that such a restric- 
tion clearly excludes the possibility of light rays traveling 
in regions denser than average. In principle, it should be 
very interesting to see how the presence of cosmic voids - 
a key entity nowadays - could be considered in the above 
prescription. 

More recently, many studies concerning the light prop- 
agation and its effects on the derived distances have been 
performed [1, 0, [3 Current constraints on the 

smoothness parameter are still weak [ll, [13, 0, Ell , how- 
ever, it is intriguing that the quoted analyses had their 
best fits for a equal to unity which corresponds to a per- 
fectly ACDM homogeneous model at all scales [l^ [l3] . 
More recently, some authors have also argued for a crucial 
deficiency of the DR approach, and, as such, it should be 
replaced by a more detailed description, probably, based 
on the weak lensing approach . 
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In this letter we advocate a slightly different but com- 
plementary point of view. It will be assumed that the 
DR approach is an useful tool in the sense that it pro- 
vides the simplest one-parametric description of the ef- 
fects caused by local inhomogeneities, but its initial con- 
ception needs to be somewhat extended. This is done 
in two steps: (i) by allowing a (here denoted by as) to 
be greater than unity in the statistical data analyses, (ii) 
by interpreting the obtained results in terms of the exis- 
tence of an uncompensated distribution of cosmic voids 
or "black regions" in the Universe (see section V). As 
we shall see, by performing a statistical analysis involv- 
ing 557 SNe la from the Union2 compilation data [HI, we 
obtain as = l-26l[J (Icr) for a flat ACDM model. This 
Icr confidence region shows that a > 1 has a very signifi- 
cant probability. We also show that a greater than unity 
is also able to harmonize the low redshift (Supernovae la) 
and baryon acoustic oscillations (BAO) data with the ob- 
servations from cosmic microwave background (CMB). 



II. THE DYER-ROEDER DISTANCE 

The differential equation driving the light propagation 
in curved spacetimes is the Sachs optical equation 



0, 



(2) 



where a prime denotes differentiation with respect to the 
affine parameter A, ^ is the cross-sectional area of the 
light beam, i?^!^ the Ricci tensor, A:^ the photon four- 
momentum {k^kf^ — 0), and the shear was neglected 20]. 

Five steps are needed to achieve the luminosity dis- 
tance in the Dyer-Roeder approach: 

• the assumption that the angular diameter distance 

• the relation between the Ricci tensor and the 
energy-momentum tensor T^i, through the Ein- 
tein's field equations 



(3) 



where in our units c = 1, i? is the scalar curvature, 
Qfii, is the metric describing a FRW geometry, G is 
Newton's constant and Rf^^kt^k" = SnGTf.^kt^k" . 

• the relation between the affine parameter A and the 
redshift z 



(4) 



where H{z) is the Hubble parameter whose present 
day value, Hq, is the Hubble's constant, 

• the ansatz pm goes to apm, and, finally. 



• the validity of the duality relation between the an- 
gular diameter and luminosity distances [13, [H, [2^ 



dL{z) - {l + zfdA{z). 



(5) 



For a general XCDM model, where the dark energy 
component is described by a perfect fiuid with equation of 
state px = wpx {w constant), the Dyer-Roeder distance 
(di = Hq^Dl) can be written as: 



aB(z)a„(l + zf + ^x{l + w){l + z)'^(i+'")] Dl{z) 

d Dl{z) ' 



-(1 



z^Eiz)^ 
dz 



(l + zYEiz) 



dz{l + z)^_ 



0,(6) 



where aE^z) denotes the extended Dyer-Roeder param- 
eter, fix, w, are the density and equation of state pa- 
rameters of dark energy while the dimensionless Hubble 
parameter, E{z) = H/Hq, reads: 

E{z) = ^n„,{l + z)3 + flxil + 0)3(1+-) + ilkil + z)2, 

(7) 

where flk — (1 — — ^x) and the limiting case 
(w = —1, ilx — ^a) of all the above expressions describe 
an arbitrary ACDM model. The above Eq.([6|) must be 
solved with two initial conditions, namely: Dl [z = 0) — 
and ^^j^\z=o = 1- As in the original DR approach, 
from now on it will be assumed that a^; is a constant 
parameter (see, however, [isl 121}). 



III. DETERMINING aE FROM SUPERNOVA 
DATA 

In order to show the physical interest of the approach 
proposed here we have performed a statistical analysis 
involving 557 SNe la from the Union2 compilation data 
[2^ . Following standard lines, we have applied the max- 
imum likelihood estimator [we refer the reader to Ref. 
[igI [2^ for details on statistical analysis involving Su- 
pernovae data]. 

In Fig. [TJa) we display the results obtained by as- 
suming a fiat ACDM model. The contours correspond 
to 68.3% (Ict) and 95.4% {2a) confidence levels. The 
best fits are ri„j = 0.25 and as — 1.26. As we can see 
from Figs. [T](b) and [TJc) the matter density parame- 
ter is well constrained, being restricted over the interval 
0.21 < ilm < 0.29 (Icr), while the smoothness parameter 
is in the interval 0.72 < as < 1.94 (Icr). Although be- 
ing poorly constrained, we see that the probability peaks 
m aE > 1, and, therefore, denser than average regions 
in the line of sight are fully compatible with the data. It 
is interesting to compare the bounds over Vim with our 
previous analysis with the restriction a < 1.0 [l6|. The 
interval 0.24 < < 0.35 {2a) was obtained. As should 
be expected, by dropping the restriction a < 1.0 lesser 
values of flm are allowed by data. 
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FIG. 1: (color online), a) The (57m, ae) plane for a flat ACDM model. The contours represent the 68.3% and 95.4% 
confidence levels. The best fit is Sim ~ 0.25 and q_b ~ 1.26. Note that a fiat model with only matter and inhomogeneities 
{Q,m = 1) is ruled out with great statistical confidence, b) Likelihood of ctB. The smoothness parameter is restricted on the 
interval 0.72 < a_B < 1-94 (Icr). c) Likelihood of fim. We see that the density parameter fi™ is restricted to the interval 
0.21 < n,„ < 0.29 (Icr). 



IV. SUPERNOVAE-CMB TENSION AND as 

The tension between low and high redshift data has 
been reported by many authors (see, for instance, [23|). 
A numerical weak lensing approach to solve this problem 
was recently discussed by Amendola et aL Q based on a 
meatball model. Can such a tension be alleviated by our 
extended DR approach? 

In order to answer that, let us consider an arbitrary 
ACDM model and plot the bounds on the {il,n,^A) 
plane by fixing three different values of a^;. By select- 
ing aE = 0.7, 1.0 and 1.3 we may study what happens 
with the (flm, ^a) contours when higher values are con- 
sidered. In Fig. we show the contours obtained for 
the chosen values of a^;. Note that when ue grows from 
0.7 to 1.3 the best fit moves of around Itr towards lower 
values of the pair (r2,„, thereby becoming more com- 
patible with the cosmic concordance flat ACDM model. 
This is a remarkable result since it improves the agree- 
ment with independent constraints coming from baryon 
acoustic oscillations (BAO) and the angular power spec- 
trum of the cosmic microwave background (CMB), and, 
more important, maintaining the same reduced Xred- 

TABLE I: Best fits for and ^a. 



as 






Xred 


0.7 


0.39 


0.83 


0.978 


1.0 


0.30 


0.78 


0.977 


1.3 


0.24 


0.74 


0.977 



In Table HI the basic results are summarized. Note 
that the greatest value of yields the minimum reduced 
Xred = X^/^ {ly is number of d.o.f). 

In Fig. EJb), we display the statistical results for a flat 
XCDM model and the same values for a e adopted in the 



previous ACDM analysis. Again, we see that for higher 
values of , the contours are displaced towards regions 
with higher values for w and smaller values for flm , again 
contributing to cancel the tension between the low and 
high redshift data. 

In Table II, we summarize the best fits for Vim and w 
along with their respective minimum reduced Xred- 

TABLE II: Best fits for SI™ and w. 



Q-E 




w 


Xred 


0.7 


0.35 


-1.18 


0.978 


1.0 


0.29 


-1.06 


0.978 


1.3 


0.23 


-0.96 


0.977 



V. WHY IS aE BIGGER THAN UNITY? 

Here we propose a simple toy model based on the exis- 
tence of cosmic voids in order to explain why ue can be 
bigger than unity. Recent studies have pointed out that 
cosmic voids not only represent a key constituent of the 
cosmic mass distribution, but, potentially, may become 
one of the cleanest probes to constrain cosmological pa- 
rameters [131 • The idea is to consider that very large 
voids are relatively rare entities, i.e. their formation suf- 
fer from the same kind of size (mass) segregation 'felt' by 
the largest galaxies and clusters. By assuming that the 3 
basic entities filling the observed Universe are: (i) matter 
homogeneously distributed (ph), (ii) the clustered com- 
ponent {pci) and (iii) voids [pvd) of small and moderate 
sizes, we define the extended DR parameter (see Eq.(l)): 



Ph + Pel + Pvd ' 
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FIG. 2: (color online) a) The influence of the smoothness parameter on the {Q.m, ^a) plane. The contours for three values of 
the smoothness parameter as using 557 SNe la from the Union2 compilation data [22| correspond to 1, 2 and 3a. Greater 
values of as provide results more compatible with a flat model, b) Contours for the {Qm , w) plane in a flat XCDM model. 
The same trend is observed, greater values of ob imply greater values of w thereby alleviating the tension among the low and 
high redshift data. 



The important task now is to quantify the contribution 
of voids representing the local underdensities in the Uni- 
verse. The presence of a void means that its matter was 
somehow redistributed to the clustered and the homo- 
geneous components. The gravitational effect of a void 
in an initially homogeneous distribution is equivalent to 
superimpose a negative density (for small densities the 
nonrelativistic superposition principle is approximately 
valid). For simplicity, it will be assumed here that the 
overall contribution of the void component can be ap- 
proximated by the linear expression, p^d — —5{ph + Pel), 
where ^ is a positive number smaller than unity. There- 
fore, ue given Eq. ([5]) can be rewritten as: 



[Ph + Pci){l - 5) 1-6 

which clearly satisfies the inequality as > a, where 
a is the standard DR parameter. In particular, when 
the clustered component does not contribute we find 
a_E = > 1. The previous analyses using supernovae 
data implies that we have effectively constrained the ex- 
tended parameter, as- How to roughly estimate the void 
contribution from this crude model? By applying the 
standard DR approach to the Union2 sample, the best 
fit is a = 1, and combining with the result for a flat 
ACDM model (section III) , one may check that the void 
contribution has a best fit of (5 ^ 0.2. It should be im- 
portant to search for a possible connection between the 
present approach and more sophisticated methods from 
weak lensing. 



VI. CONCLUSIONS 

In this Letter we have discussed the role played by lo- 
cal inhomogeneities on the light propagation based on 
an extended Dyer-Roeder approach. In the new inter- 
pretation light can travel in regions denser than average, 
a possibility phenomenologically described by a smooth- 
ness parameter > 1. 

In order to test such a hypothesis we have performed 
a statistical analysis in a flat ACDM model and the best 
fit achieved was as = 1.26 and flm = 0.25, the param- 
eters being restricted to the intervals 0.72 < < 1.94 
and 0.21 < < 0.29 within the 68.3% confidence level. 
Although as being poorly constrained, the results are 
fully compatible with the hypothesis of light traveling 
in denser than average regions. We have also analyzed 
how different values for the smoothness parameter affect 
the bounds over {flm, ^a) in an arbitrary ACDM model. 
Interestingly, > 1 improves the cosmic concordance 
model since it provides a better agreement between low 
and high redshift data (Supernovae, CMB and BAO). 
The same happened when a flat XCDM model was con- 
sidered with the assumption that > 1. 

Such results suggest that the hypothesis of light travel- 
ing in regions denser than the cosmic average seems to be 
quite realistic. A toy model justifying why this may occur 
with values of greater than unity was also discussed 
by taking into account the possible influence of cosmic 
voids on the Dyer-Roeder approach. The simplicity of 
the model and the obtained results reinforce the interest 
on the influence of local inhomogeneities and may pave 
the way for a more fundamental description. 
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